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ABSTRACT 

In the standard cosmological model, the temperature anisotropy of the cosmic 
microwave background is interpreted as variation in the gravitational potential 
at the point of emission, due to the emitter being embedded in a region C of over- 
or under-density spanning the length (or size) scale A on which the anisotropy is 
measured. If the Universe is inhomogeneous, however, similar density contrasts 
of size A are also located everywhere surrounding C. Since they are superposition 
states of many independent Fourier modes with no preferred direction, such pri- 
mordial clumps and voids should not be configured according to some prescribed 
spatial pattern. Rather, they can randomly trade spaces with each other while 
preserving the Harrison- Zeldovich character of the matter spectrum. The out- 
come is an extra perturbation of the potential when averaged over length A at the 
emitter, and consequently an additional anisotropy on the same scale, which has 
apparently been overlooked. Unlike the conventional application of the Sachs- 
Wolfe effect to the WMAP observations, this extra effect is not scale independent 
over the P(k) ~ k part of the matter spectrum, but increases towards smaller 
lengths, as yk. The consequence is a substantial revision of the currently ad- 
vertised values of the key cosmological parameters, unless one postulates a more 
rapid decrease in the gravitational force with distance than that given by the 
inverse- square law. 

1. Introduction and preliminaries 

In their seminal paper, Sachs and Wolfe (1967) discussed the effect of an inhomogeneous 
Universe at the decoupling epoch on the temperature anisotropy of the cosmic microwave 
background (CMB). They considered random density perturbations Sp/p over some charac- 
teristic length scale A - a phenomenon which was subsequently conjectured to have originated 
from scale invariant fluctuations in the primordial plasma (Harrison 1970, Zeldovich 1972, 
Peebles & Yu 1970). This behavior, and its manifestation as CMB temperature anisotropy, 
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has been investigated in detail both theoretically (Peebles 1982, Bond & Efstathiou 1984 and 
1987) and observationally (Bennett et al 2003, Page et al 2003, and Spergel et al 2003 pro- 
vided results from the latest all-sky measurements by WMAP). The purpose of the present 
work is to point out that the conversion from matter inhomogeneity to CMB anisotropy, 
adopted by conventional interpretation of the WMAP TT cross correlation data, is a re- 
stricted version of the Sachs- Wolfe effect which does not take into account all contributions 
(to the anisotropy) of comparable magnitude. 

Let us first set up the necessary preliminary framework. At a spatial corrdinate r and 
referring all physical quantities to their values at the present epoch where the expansion 
parameter is a = 1 let the matter density be p(r) = p + 5p(r) with the perturbation term 
having zero spatial average over some large volume, i.e. (5p(r)) = 0. The Fourier transform 
of p, and its inverse transform, are 

p(k) = I dhe-^p(r), p(v) = [ 7^e^p(k), (1) 



(27T)* 

respectively. The deviation in the gravitational potential at r, due to the inhomogeneous 
Universe, is given by 

6*(r) = -Gj dh'-^l. (2) 

We shall return to Eq. (2) soon. For the moment, let us examine 5p(k), because the power 
spectrum P(k) is characterized by an equation which involves 5p: 

MVW) = P lP(k)(2nf5(k + k'). (3) 

For the Harrison- Zeldovich (HZ) spectrum, P(k) ~ k. 

We now define a smoothing function (or filter) W\(r) corresponding to any chosen scale 
A, satisfying W\(r) ~ 1 for r <C A and W\(r) ps for r ^> A. The effective smoothing volume 
V x is: 

V x = J dhW x {v). (4) 

For example, a convenient choice is W\(r) = e ~ r2 / 2A2 for which V\ = (27r) 3 / 2 A 3 and W\(k) = 
V x e~ x k I 2 . Then the smoothed mass over length scale A will be 

5M x (r) = J dh' W\{r - r')5p(r'), (5) 

or equivalently 

5M x (k) = W x (k)5p(k)- (6) 



- 3- 



The averages of M x are (5M x {v)) = and 

j^ 3 e^-^\W x (k)fplP(k). (7) 

Likewise for the smoothed density Sp x we have (5p x (r)) = and 

(5~p x (k)5p x (V)) = pte- X2k2 P(k)(2nf5(k + k'). 

In particular when P(k) = Ak n , the mean square value of the smoothed mass can be eval- 
uated by the saddle-point approximation. The k integral is dominated by the region near 
k — k\— y/1 + f/A, to give 

{{5M x f) = c n M 2 x klP(k x ), (8) 

where M x = PqV x and c n is a numerical constant. Here-and-after we shall define the standard 
deviation in 5M x (r) as 

5Mr = ((5M x ) 2 )K (9) 

A note of caution is already in order here. Although the r.m.s. mass and density over some 
length scale A at position r concern the distribution of matter local to r, the same cannot be 
said about r.m.s. values of the potential fluctuation 5<&(r), Eq. (2). Owing to the long range 
nature of the gravitational force, density contrasts of size A but spreading over distances far 
greater than A can also contribute towards <5$™ s . 

2. CMB temperature anisotropy from primordial density fluctuations: any miss- 
ing component? 

If at some point r on the last scattering surface there is a mass excess (say) of <5M™ S 
over length scale A, the most obvious contribution to the CMB anisotropy in this scale will 
be a perturbation in the gravitational potential of the form 

G5M X G5M X 
5$ A (r) w w -^- Po V x . (10) 



Since V x ~ A 3 and from Eq. (8) we have 5M X /M X ~ k\ ~ A 2 for P(k) ~ k, it is then 
apparent that 

~ <^ = con stant, (11) 

i.e. for primordial HZ fluctuations the value of smoothed temperature anisotropy is indepen- 
dent of the length scale of the smoothing filter. This is the well known Sachs- Wolfe effect as 
applied to a HZ matter spectrum. 



Let us however query whether the standard result described above, which concerns 
anisotropics caused by fluctuation in the mass surrounding the site of photon emission and 



-4- 



extending to the length scale A under consideration, represent all the possibilities, Of course, 
on scales A any total mass variation will, in the restricted context of Eqs. (10) and (11), 
lead to anisotropies over correspondingly larger angular separations. Yet Eq. (10) is not the 
only way to perturb $A( r )- A simple analogy is that the potential at some point on the earth 
surface depends solely on the total mass of the earth if there is perfect spherical symmetry in 
the density function. Deviations in $ from point to point on the surface can occur even when 
the total mass is fixed, if non-uniformities in the matter distribution exists in any part of 
the earth's interior, far or near the location in question. Moreover, the distance over which 
potential excursions occur equals the mean spacing between mass concentrations. In the 
present problem, obviously beyond the distance A from a CMB emitter at position r, space 
remains just as inhomogeneous on the scale of A. Then, different configurational realizations 
(or random placement) of these size A density contrasts (here-and-after referred to as lumps) 
within the, horizon at decoupling, i.e. the causal sphere TZ of radius R ^> A centered at the 
emission point r, can also result in a finite <5$a(]~) which is unrelated to Eq. (10). 

Thus our contention is that when the distribution of primordial matter is smoothed at 
resolution A, the resulting lumps of over- and underdensity, being superposition states of 
many independent Fourier modes with no preferred direction, can randomly trade spaces 
with each other and the HZ nature of the matter spectrum (on scales upwards of A, of 
course) will still be preserved 1 . This leads to an additional component of 5$>(r) from 
the configurational arbitrariness of an entire ensemble of lumps, all in causal contact with 
position r. 

It is in fact quite easy to estimate the magnitude of the additional contribution. If 
primordial lumps of various masses and size A randomly pack the Universe with 100 % filling 
factor, the resulting total mass fluctuation in any multiple-lump sub-region will remain in 
compliance with 5M/M ~ k 2 , yet there will be a perturbation on the potential &\(r), given 
by 5$™ s /$a ~ l/y/N where N as i? 3 /A 3 is the total number of lumps in our causal sphere, 
although we shall derive the precise value of <5$™ s in the next section. Thus we deduce that 
there is more CMB temperature anisotropy to be reckoned with, of order 



which has not been taken into account by conventional (single lump) application of the 
Sachs- Wolfe effect, Eq. (11). For a HZ matter spectrum, use of Eq. (8) enables us to write 



If the matter is of primordial origin the fluctuation in the total mass of any multiple lump region will 
remain consistent with the HZ spectrum, irrespective of how the lumps are re-arranged, provided that the 
volume filling factor (by the lumps) is 100 %, and the procedure of re-configuration ensures randomness of 
the lump positions within some very large volume. 




(12) 
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Eq. (12) as 

'ST\ IR 



- J ~ constant \ — , (13) 
/a V A 

where the constant factor is the same as that in Eq. (11). Bearing in mind that R > A, the 
new effect presented here appears at least as important as the standard result. 

3. Anisotropy from analysis in configurational space 

A slightly different way of looking at the physics described in section 2 is afforded by 
working in real space. After spatial filtering the excursions in the potential, Eq. (2), can be 
written as a discrete sum over the resolution elements, viz. 

«*)=E^nr (14) 

Here the mass fluctuation at position i is given by SM{, a quantity with zero mean and 
standard deviation SM r x ms where 

SM r x ms = ^7rA 3 5p r A ms . (15) 

Now the contribution to 5<l>,\(r) from the density contrast in the vicinity of the point r, i.e. 
the 'local lump' at r ; = r where |r — r | ~ A, is 

„ , , G5M T x ms G 5M r x ms 4 ^ 1 TT2 , 2 5p x ms 

where Q m is the mean matter density of the present Universe in units of the critical den- 
sity p c = 3Hq/(8ttG). This yields a temperature change (as the emitted CMB radiation 
leaves its own region of potential excursion) of ST/T as 5<&\, in agreement with Eq. (48) 
of Sachs & Wolfe (1967). Note that because in the case of a HZ spectrum the quantity 
Sp r x ms /p\ ~ <5M™ s /Ma ~ 1/A 2 , we have from Eq. (8) 5<&\(r) = constant, in agreement 
with the conclusion of the previous section that single lump Sachs- Wolfe effect leads to the 
standard result of scale invariant CMB anisotropy for primoridal matter distributions. 

The important point, however, is that Eq. (14) also depicts effects beyond that of the 
local lump. As stated earlier, the contribution to <5$A(r) from the random placement all 
other lumps within the horizon should be included with the summation procedure. It is 
emphasized again that this effect is not because of the Poisson statistics in the total mass 
of many independently varying density contrasts belonging to the causal sphere 1Z (which 
is suppressed by the HZ spectrum of P(k) ~ k), but because of the Poisson process in the 
configurational arrangement of the same primordial lumps inside 71. 
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Our task initially is to compute the deviation <5$™ s in the gravitational potential of a 
sphere of radius R when it contains a constant number N of smaller spherical lump of radius 
A, each randomly placed and filled with matter accounting for a mass of 5M" ns per lump, 
such that the remaining matter-free regions occupy half the volume of the causal sphere. 
Next, it is realized that this actually means neglecting the contribution to 5$" ns from the 
underdense regions (each of mass — <5M™ S ), when it is included <5$™ s will increase by a factor 
of two because the gravitational effects of these two types of regions are anti-correlated. The 
potential deviation at the position r (r < R) due to the overdense lumps is 

*Mr) = -£Tn TP ( 17 ) 

i 

Now each of the N overdense lumps has the same probability distribution, in which the 
probability p(r)d 3 r of finding the lump in a small volume <i 3 r is given by 

p(r)dh = ^dh. (18) 

where in Eq. (18) we adopted the continuum approximation appropriate to the limit of many 
lumps. Hence the mean potential ($(r)) is 

(6$ x (r)) = -NG5MI™ ( = -GSMT 8 [ n ^ d ^' (19) 

\ |r — r I / J |r — r I 

Relevant to the present problem is a fixed set of overdense lumps within TZ, i.e. n(r) = 
37V/ (47T.R 3 ) for r < R, and n(r) = for r > R since the matter lying beyond 1Z plays no 
part. The integral of Eq. (19) can readily be evaluated under this scenario to yield 

where Eq. (20) is valid for the range < r < R. Thus at r = it! one recovers the usual New- 
tonian potential (6&\) = —GN5Ml ms /R, which on average cancels exactly the contribution 
from the underdense lumps. 

To obtain the variance {(5Q\) 2 ), the quantity by which the potential for points in space 
separated by lengths > 2A differ from each other, we need the mean square of <5$a defined 

(W')]VE(S5). PI) 

When we subtract from this the 'square of the mean', viz. the quantity 
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it is clear that all the terms with j ^ k will cancel, and we are left with 



2 X 




<[M>A(r)] 2 > - (5$A(r)) 2 = G\5MD 2 N — - — — . (23) 



Now the average related to the first term on the right side of Eq. (23) may also be computed 
analytically for the form of n(r) already discussed. By this we mean 



/ 1 \ 3 r r'dr' ( r + r' \ 
\|r-r'| 2 / ~ 2R 3 J Q 11 \\ r -r'\) 



(24) 



At the surface itself r = R, and the integral reduces to the simple form 

1 \ 3 



at r = R. (25) 



- - r'| 2 / 2R 2 

The 2nd term on the right side of Eq. (23) is, by Eq. (19), equal to (5$ X ) 2 /N. 

Thus, altogether and after including the aforementioned contribution to 5$™ 1S from the 
underdense lumps, we arrive at a standard deviation of 



5$ r A ms (r) 



3G 2 (5M[ ms ) 2 N G 2 (5M[ ms ) 2 N 
2R 2 R 2 



f 2N G5M{ ms 



R 



(26) 



Eq. (26) depicts the spatial variation in the potential between two points located at a 
distance > 2A apart. 

The CMB anisotropy is finally obtained from Eq. (26) by observing that within the 
causal sphere there are altogether R 3 /\ 3 lumps of radius A, half of which are overdense. 
Thus N = R 3 /{2\ 3 ), and we have 



rms 



^5§ T ™ = G5Ml ms \-lR^ = -n m H 2 \lR^^- for X < R. (27) 
T J x ' 2 Pa 

This gives the new contribution which should explicitly be invoked 2 as a separate anisotropy 
term in the 'Boltzmann solver' codes like CMBFAST (the routine used to fit the standard 
cosmological model to WMAP data, see Zaldarriaga & Seljak 2000), but is not. More 
elaborately, for self-consistency of the standard model this contribution must be included, 
even though it has hitherto been ignored. When compared with Eq. (16), we see that 
although conventional Sachs- Wolfe effect yields a constant 5T/T at all scales, the extra 



2 That is to say, the term will not arise 'naturally' as solution of a large number of coupled differential 
equations unless the physics behind it is in place. 
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anisotropy (5T/T)\ equals this same constant multiplied by the factor y/R/X (in agreement 
with the earlier analysis of section 3), i.e. it is not scale invariant. The reason why CMB 
temperatures must exhibit the variation of Eq. (27) is that two points on the last scattering 
surface separated by a distance > 2A experience statistically independent potential deviation, 
due to the different realization (or pattern) of lumps in the causal spheres centered at these 
points. 

It should also be mentioned that no matter how many more contrainsts unassumed 
by us were to be applied to restrict the number of permitted lump configurations within 
7Z, so long as space is not totally homogeneous beyond the emitter's own 'local lump' the 
incorporation of effects of remote density contrasts via the long range force of gravity must 
remove the constancy of (8T/T)\ by introducing a factor that scales with R/\ in some way. 
The fact that the standard model predicts a constant (5T/T)\ at large (but sub-horizon) A 
offers the clearest indication of it's failure in recognizing the phenomenon presented here. 

4. Have the CMB anisotropy observations been interpreted correctly? 

At sufficiently low spherical harmonics i < 90 where the smoothing length A > R, lack of 
causal linkage between lumps renders Eq. (27) inapplicable, and the only anisotropy would 
be the conventional single lump Sach- Wolfe effect of Eq. (16). From the WMAP TT cross 
correlation plot (Bennett et al 2003) one sees that (5T/T)\ m 10~ 5 at £ <C 90. Thus we 
conclude that 

^ m tf 2 A 2 ^™ » 10- 5 (28) 
^ Pa 

on all scales A so long as the matter distribution follows the primordial HZ spectrum of 
P(k) ~ k. 

For sub-horizon features at decoupling which were super-horizon in size during matter- 
radiation equipartition 3 , like the first acoustic peak, the extra anisotropy of Eq. (27) must 
now be included. Since R ~ 490 Mpc (present value for the causal radius at decoupling) and, 
in the case of the first peak the full size of these lumps is 2A ~ 147 Mpc, we see from Eq. 
(27) and (28) that the total expected anisotropy with the conventional and new component 
added in quadrature is 

ST 

— « 2.77 x 10~ 5 (29) 



3 Structures that fit inside the comoving equipartition horizon have according to the standard model 
a density contrast commensurate with P(k) ~ k~ 3 at the primordial level. Thus our present conclusion 
regarding the first peak, which is however based upon a P(k) ~ k spectrum, does not apply to the higher 
peaks. 
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for the first peak. Obviously, Eq. (29) in conjunction with the WMAP observations pose a 
dilemma for the standard cosmological model, because the detected first acoustic anisotropy 
is only pa 2.5 x 1CT 5 , it already equals the expected value as given by Eq. (29), which is based 
upon a purely primordial spectrum before enhancement of density contrasts by sound waves. 
This leaves no room for interpreting the acoustic features as compression and rarefaction of 
the plasma by sound propagation. 

Is there a way of restoring the standard cosmological model, which matches the WMAP 
data very well? We could supposedly alter the cosmological parameters, e.g. the primordial 
spectral index could substantially increase from n — 1 to reduce the impact of Eq. (27), 
though the epistimological justification of such changes will almost certainly be artificial. 
One less damaging possibility remains, however. The gravitational force beyond some > 100 
Mpc distance scale could wane more quickly than the inverse-square law. This would appear 
to the author as a remedy of 'least resistance' for several reasons. Firstly, the standard 
model can immediately be reinstated, as the Friedmann equations are constructed in the 
context of General Relativity without appealing to the existence of a gravitational influence 
on Hubble scales. Secondly, structure formation theory is unaffected, because numerical 
codes that simulate the building galaxies, groups, and even clusters do not depend on the 
inverse- square law persisting to beyond 100 Mpc. Thirdly, there is currently no reliable 
experimental constraint on the behavior of gravity at these very large distances. 
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